An Analytical Solution
Federal Reserve Bank of Kansas City
KC Fed Brown Bag
August 4 2026
The views expressed herein are those of the author and do not necessarily represent the views of the Federal Reserve Bank of Kansas City or the Federal Reserve System.
Macroeconomics with fixed costs of adjustment depends on the distribution of agents in the inaction region
This requires solving a PDE: the Kolmogorov Forward Equation (KFE)
Analytical solutions aid tractability and understanding (e.g. Alvarez–Lippi sufficient statistics)
But an analytical solution is difficult because the PDE is endogenous: its evolution depends on the flow of resets
Existing methods require shortcuts, such as symmetry or small shocks
Derive an analytical solution to a generic fixed cost model (KFE + boundary conditions)
Key insight: in the frequency domain, you can solve for the endogenous flow of resets without first finding the entire distribution
Reset flow is like magic! All aggregate IRFs are linear functions of the reset flow alone: entire distribution is not needed
Opens up new theoretical results!
Application: the advent of AI with investment adjustment costs
Examples: sticky prices, investment, inventories, rational inattention, hiring and firing, wage negotiation, etc.
Agent’s state \(x\) follows a diffusion inside inaction region; agent pays a cost to adjust \(x\) only occasionally
\[ dx=\mu\,dt+\sigma\,dW, \qquad x\in[0,\bar x], \qquad x\longrightarrow x^* \]
Agent’s only decision is the inaction region + optimal reset. Well-understood problem!
Aggregation (macroeconomics) is the challenge
The density \(h(x,t)\) follows the Kolmogorov Forward Equation (KFE): \[ \partial_t h(x,t) =\frac{\sigma^2}{2}\,\partial_x^2 h(x,t)-\mu\,\partial_x h(x,t)-\eta h(x,t) +\delta(x-x^*)F(t) \] \[ h(0,t)=h(\bar x,t)=0 \]
Lemma 1: flow \(F(t)\) of agents resetting at time \(t\) is
\[ F(t) =\frac{\sigma^2}{2}\,\partial_x h(0,t) -\frac{\sigma^2}{2}\,\partial_x h(\bar x,t) +\eta \]
The known Green’s function \(G(x,y,t)\) is the response of density over \(x\) at \(t\) from impulse at \(y\)
Textbook formula for \(h(x,t)\) conditional on \(F(t)\):
\[ h(x,t) = \underbrace{\int_0^{\bar x}h(y,0)G(x,y,t)\,dy}_{\text{initial distribution}} + \underbrace{\int_0^tF(\tau)G(x,x^*,t-\tau)\,d\tau}_{\text{earlier reset cohorts}} \]
Theorem 1: Expression for \(\widehat{F}(s)\) (Laplace transform) in terms of parameters, initial condition: \[ \widehat F(s)=\frac{\widehat A(s)}{1-\widehat B(s)} \] \[ F(t)=A(t)+\int_0^t B(t-\tau)F(\tau)\,d\tau \]
\[ Z(t)=\int_0^{\bar x}f_Z(x)h(x,t)\,dx \]
\[ \widehat Z(s)=\widehat A^Z(s)+\widehat B^Z(s)\widehat F(s) \]
\[ Z(t)=A^Z(t)+\int_0^tB^Z(t-\tau)F(\tau)\,d\tau \]
Closed form solutions apply to the (common) baseline
But insights (importance/utility of reset flow) are more general
\[ \partial_t h(x,t) = \partial_x^2\!\left(\frac{\sigma^2(x)}{2}h(x,t)\right) -\partial_x\!\left(\mu(x)h(x,t)\right) -\eta(x)h(x,t) +\delta(x-x^*)F(t) \] \[ h(0,t)=h(\bar x,t)=0 \]
\[ F(t) = - \left.\partial_x\!\left(\frac{\sigma^2(x)}{2}h(x,t)\right)\right|_{x=0}^{x=\bar x} + \int_0^{\bar x}\eta(x)h(x,t)\,dx \]
Analytical solution opens possibility of useful new theoretical results
Example: a new sufficient statistic approach for the effects of macro shocks
Existing results map steady-state adjustment moments into the cumulative impulse response (e.g. Alvarez-Lippi-Le Bihan 2016)
Specific to some specialized settings, particular aggregates (e.g. the avg. state)
Can we learn more?
\[ S(t\mid y)=\int_0^{\bar x}G(x,y,t)\,dx \qquad m(t\mid y)=-\partial_tS(t\mid y) \] \[ \acute{\mathbf S}(t) \equiv\int_0^{\bar x}\partial_yS(t\mid y)\,\bar h(y)\,dy \qquad \mathbf m(t)\equiv m(t\mid x^*) \]
\[ \widehat F^\Delta(s) = \frac{-s\,\widehat{\acute{\mathbf S}}(s)} {1-\widehat{\mathbf m}(s)} \]
\[ F^\Delta(t) =-\partial_t\acute{\mathbf S}(t) +\int_0^t\mathbf m(t-\tau)F^\Delta(\tau)\,d\tau \]
Steady state objects \(\mathbf m\) and \(\acute{\mathbf S}\) identify the entire marginal response path
\(F^\Delta(t)\) + Theorem 2 \(\longrightarrow\) get the marginal IRF of other aggregates
A special case of Baley and Blanco (2021); KFE is:
\[ \partial_t h(x,t) = \frac{\sigma_p^2}{2}\,\partial_x^2h(x,t) +(\zeta+g_p)\,\partial_xh(x,t) +\delta(x-x^*)F(t) \qquad h(0,t)=h(\bar x,t)=0 \]
Common to focus on average gaps as approximations to other aggregates
Why? \(h(x,t)\) gives dynamics for avg \(\tfrac{K}{P}\)
But joint dist. \(f(p,x,t)\) is needed for \(K(t)\), \(Y(t)\), etc.
… or is it? Another insight from analytical solution: we do not need the joint distribution to derive aggregate IRFs!
\[ \begin{aligned} K_P(t)\equiv\frac{K(t)}{\bar P(t)} &=\int_0^{\bar x}e^{\underline{\widetilde k}+x}h_P(x,t)\,dx\\ Y_P(t)\equiv\frac{Y(t)}{\bar P(t)} &=\int_0^{\bar x}e^{\alpha(\underline{\widetilde k}+x)}h_P(x,t)\,dx \end{aligned} \]
\[ \partial_t h_P =\frac{\sigma_p^2}{2}\partial_x^2h_P + \left(\zeta+g_p + \sigma^2\right)\partial_xh_P +\delta(x-x^*)F_P(t) \]
\[ F_P(t)=\frac{\sigma_p^2}{2} \left(\partial_xh_P(0,t)-\partial_xh_P(\bar x,t)\right) \]
\(h_P(x,t)\) inherits absorbing boundaries, integrates to one, all density re-enters at \(x^*\)
The reset-flow solution therefore applies to \(F_P(t)\) and \(h_P(x,t)\)
Can also use Theorem 2 to derive aggregates, Theorem 3 to calculate the suff. stat.
Analytical shortcut: can solve for the endogenous reset flow on its own (Theorem 1), which gives you everything else
Aggregation: many macroeconomic variables are linear functions of the reset flow (Theorem 2)
Empirical implication: conditional adjustment behavior identifies the full marginal reset-flow IRF (Theorem 3)
Investment: a productivity-weighted gap density is sufficient for aggregate capital, output, and investment
Economic result: faster productivity growth generates an immediate investment boom followed by an extended lull whose depth depends nonlinearly on shock size
Derivations, validation, and robustness
\[ \partial_t h =\partial_x^2\!\left(\frac{\sigma^2(x)}{2}h\right) -\partial_x\!\left(\mu(x)h\right) -\eta(x)h +\delta(x-x^*)F(t) \]
\[ F(t) = \left.\partial_x\!\left(\frac{\sigma^2(x)}{2}h(x,t)\right)\right|_{x=0} - \left.\partial_x\!\left(\frac{\sigma^2(x)}{2}h(x,t)\right)\right|_{x=\bar x} + \int_0^{\bar x}\eta(x)h(x,t)\,dx \]
\[ G(x,y,t) =\sum_{n=1}^{\infty} \psi_n^r(x)\psi_n^\ell(y)e^{-\lambda_nt} \]
\[ \widehat A(s) =\sum_{n=1}^{\infty}\frac{\theta_na_n}{s+\lambda_n}, \qquad \widehat B(s) =\sum_{n=1}^{\infty}\frac{\theta_nb_n}{s+\lambda_n} \]
\[ a_n=\int_0^{\bar x}\psi_n^\ell(y)\phi(y)\,dy, \qquad b_n=\psi_n^\ell(x^*) \]
Value matching
\[ \begin{aligned} v(\underline{\widetilde k}) &=v(\widetilde k^*)-q \left(e^{\widetilde k^*}-e^{\underline{\widetilde k}}\right)-\kappa\\ v(\bar{\widetilde k}) &=v(\widetilde k^*)-q \left(e^{\widetilde k^*}-e^{\bar{\widetilde k}}\right)-\kappa \end{aligned} \]
Smooth pasting and optimal reset
\[ \begin{aligned} v'(\underline{\widetilde k})&=q e^{\underline{\widetilde k}}\\ v'(\bar{\widetilde k})&=q e^{\bar{\widetilde k}}\\ v'(\widetilde k^*)&=q e^{\widetilde k^*} \end{aligned} \]